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ABSTRACT 



43 

Supersymmetric Unified theories which incorporate a renormalizable Type I seesaw 
^■l mechanism for small neutrino masses can also provide slow roll inflection point inflation 

along a flat direction associated with a gauge invariant combination of the Higgs, slepton 
and right handed sneutrino superfields. Inflationary parameters are related to the Major ana 
and Dirac couplings responsible for neutrino masses with the scale of inflation set by a right- 
handed neutrino mass M u c ~ 10 6 — 10 12 GeV. Tuning of the neutrino Dirac and Majorana 
superpotential couplings and soft Susy breaking parameters is required to enforce flatness of 
the inflationary potential. In contrast to previous inflection point inflation models the cubic 
term is dominantly derived from superpotential couplings rather than soft A-terms. Thus 
since M u c >> Ms usy the tuning condition is almost independent of the soft supersymmetry 
breaking parameters and therefore more stable. The required fine tuning is also less stringent 
than for Minimal SUSY Standard Model (MSSM) inflation or Dirac neutrino "A-term" 
inflation scenarios due to the much larger value of the inflaton mass. Reheating proceeds 
via 'instant preheating' which rapidly dumps all the inflaton energy into a MSSM mode 
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radiation bath giving a high reheat temperature T rh » M„ 4 C 10 6 GeV ~ 10 11 - 10 15 GeV. 
Thus our scenario requires large gravitino mass > 50 TeV to avoid a gravitino problem. 
The 'instant preheating' and Higgs component of the inflaton also imply a 'non-thermal' 
contribution to Leptogenesis due to facilitated production of right handed neutrinos during 
inflaton decay. We derive the tuning conditions for the scenario to work in the realistic New 
Minimal Supersymmetric SO(10) GUT and show that they can be satisfied by realistic fits. 
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1 Introduction 



Primordial inflation is perhaps the simplest dynamical mechanism which can explain the 
seed fluctuations pQ for the cosmic microwave background (CMB) radiation, and thus for the 
formation of large scale structures. Although a large number of inflationary models exist 
in the literature the majority of them are not grounded in any realistic model of particle 
physics, thus leaving them unconstrained by anything beyond the few parameters so far 
gleaned from measurements [21 [3] of the Cosmic Microwave background. Models where 
inflation is driven not by a generic scalar field but by an inflaton intimately tied[l] to the 
Standard Model gauge group and spectrum carry an obvious appeal. Moreover in order to 
have a successful (and calculable) reheating into the Standard Model degrees of freedom 
as is required for the success of Big Bang Nucleosynthesis, the inflaton must carry definite 
Standard Model gauge and Yukawa charges so that the inflaton condensate can efficiently 
decay into SM degrees of freedom after the end of inflation. 

The suggestion [3] that inflation can be naturally embedded within the Minimal su- 
persymmetric (SUSY) Standard Model (MSSM), with generic gravity mediated (i.e N = 1 
supergravity type : but we assume canonical Kahler potential ) soft supersymmetry break- 
ing terms, is an attractive scenario which enables us to connect the microscopic origin of 
inflation to cosmological evolution on the largest scales. Models of this type are typically 
based on slow roll inflation associated with "flat directions" in the MSSM field space (along 
which the D-term potential vanishes). A well known theorem [5] allows one to use holomor- 
phic gauge invariants formed from chiral superfields as coordinates for the D-flat manifold 
of the scalar field space of SUSY gauge theories. The flat directions are lifted by super- 
gravity generated soft supersymmetry breaking terms and by non renormalizable terms in 
the MSSM effective superpotential. Such models ( also called "A- Term Inflation" models 
[SJ El E]) typically require a fine tuning between the soft terms to ensure an inflection or 
saddle point of the field potential where the vacuum energy density drives a burst of in- 
flation but nevertheless allows "graceful exit" due to the absence of a local minimum and 
the associated potential barrier which would prevent exit. In such models the (usually non- 
renormalizable) terms that lift D-flatness of the inflaton potential are hypothesized rather 
than deduced from a well defined underlying renormalizable model. Thus while they answer 
some of the relevant issues they have much scope for improvement. One may consider how 
to deduce the effective non-renormalizable superpotential by integrating out heavy fields 
from an underlying theory, or one may look for minimal extensions of the MSSM which 
may (like inflationary GUTs) support inflation even at the renormalizable level. 

The first definite signal of physics beyond the SM came from neutrino oscillations which 
are now accepted as evidence of non zero neutrino masses in the milli-eV range. However 
the nature of neutrino masses, i.e whether they are of Dirac or Majorana type, is still 
unsettled. In the first case light neutrino masses are understood as being the consequence 
of highly suppressed Yukawa couplings, (£ = y v N rH u L + y u ~ m^/Mw ~ 10 -12 ) 7 
or more orders of magnitude smaller than the charged fermion Yukawa couplings. To be 
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dominantly of Dirac type these masses should be accompanied by highly suppressed right 
handed Majorana neutrino masses, M u c ~ O.leV or less. Conversely one may generate 
small effective (Type I seesaw[9]) neutrino Majorana masses (m v ~ (m{?) 2 /M u c) for the left 
handed neutrinos if the right handed neutrino masses M u c take the large values permitted 
by their vanishing SM gauge charges. In this case the Dirac masses of the neutrinos need 
not be suppressed by ultra small Yukawa couplings as required in the Dirac mass case. 

In [6j[7] an intriguing connection was made between the smallness of the (Dirac) neutrino 
masses and flatness of the inflaton potential within the MSSM extended by the addition of 
U(1)b-l gauge group and right handed neutrinos. The inflaton field was a gauge invariant 
D-flat direction, NH U L, where N is the right handed sneutrino, H u is the MSSM Higgs 
which gives masses to the up-type quarks, and L is the slepton field. The gauge invariant 
superpotential term y v NH u L generates the tiny (Dirac) neutrino masses due to the afore- 
mentioned tiny neutrino Yukawa coupling (y v ~ 10~ 12 ). When coupled with soft trilinear 
and bilinear supersymmetry breaking terms of mass scale ~ 100 GeV to 10 TeV the asso- 
ciated renormalizable inflaton potential can then be fine tuned to achieve inflection point 
inflation consistent with Wilkinson Microwave Anisotropy Probe(WMAP) 7 year datajfJJE]. 

Since the seesaw[9] explanation for neutrino masses is arguably preferable to the ad-hoc 
small Dirac masses explanation it is natural to ask if it too supports inflation. Prima facie 
such a scenario could face obstacles in meeting the requirements of the neutrino- inflaton sce- 
nario i.e ultra small superpotential couplings, and TeV scale trilinear /mass terms. Generic 
Type I seesaw relies upon large right handed neutrino Majorana masses which are gener- 
ated by breaking of B — L symmetry by vevs Vb~l » 10 10 Gey. An inflaton involving the 
right handed sneutrino will then have (supersymmetric) mass contributions as large as the 
righthanded neutrino mass. The cogency of the seesaw lies in not artificially singling out 
neutrino Yukawas to be ultra small. With large Vb-l the Dirac coupling of the neutrino 
need not be suppressed by hand. Indeed for, normal hierarchy, one obtains the third gen- 
eration light neutrino masses m u ~ y 2 v 2 EW / M v c ~ O.leV for 7/3 ~ 1, M u c ~ 10 15 GeV. Such 
large couplings and masses would completely destroy the needed flatness of the inflationary 
potential. However reflection shows that a negative conclusion may be unwarranted since 
at least three generations of neutrinos and their superpartners are in play. So there is con- 
siderable scope for much smaller superpotential couplings: the neutrino Yukawa coupling 
eigenvalues could have the typical values associated with up type fermions while off diagonal 
components matched the tiny Majorana couplings in smallness. Off diagonal flat directions 
(NaHLb,A 7^ 5=1,2,3) can serve just as well as diagonal ones , in fact we shall see they 
are required in the realistic New Minimal SO (10) GUT implementation of our scenario . 

Furthermore the popular Leptogenesis[10 scenario strongly hints at right handed neu- 
trino masses in the range 10 6 to 10 12 GeV. So for Vb-l ~ Mx > 10 16 GeV the superpoten- 
tial couplings Ja, A = 1,2,3 (we will work in a basis where these couplings are diagonal), 
which generate right handed neutrino masses M u <^ ~ JaVb-l, are very small (/a ~ 10~ 9 
to 10 -4 ). Thus the required ingredients for an inflaton in the Type I seesaw scenario are 
already present. Note that since generic Type I seesaw requires that B — L is broken at 
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a high scale, issues concerning the efficient decay of the conjugate sneutrino component 
of the inflaton via their Yukawa couplings will need to be addressed. Reheating in our 
scenario proceeds via the so called 'instant pre-heating' mechanism [11] resulting in a high 
reheat temperature due to rapid dump of the innaton energy into MSSM modes. The Higgs 
component of the innaton implies [12] a non-thermal contribution to leptogenesis. 

Issues regarding natural values for superpotential couplings come into focus when viewed 
in the context of the so called Minimal Left Right supersymmetric models [13j and their em- 
bedding in GUT models [T4"l [T5] . SUSY Left-Right Models are advantaged due to their 
protection of R-parity as a gauged discrete symmetry, which provides a stable lightest su- 
persymmetric particle (LSP). They simultaneously and naturally implement Seesaw mech- 
anisms for neutrino masses [IB]. Moreover such models have also been incorporated in the 
realistic and predictive New Minimal Susy SO(10) grand unified theories(NMSGUT) [16} [T7] 
where all the hard parameters of the MSSM are fitted in terms of fundamental parameters of 
the GUT and soft SUSY breaking parameters (of the Non-Universal Higgs masses (NUHM) 
type) defined at the Unification scale Mx ~ 10 16 — 10 18 GeV. Such GUTs have viable 
Bino dark matter candidates and make distinctive predictions for the type of SUSY spec- 
tra observable at the LHC. In 2008, well before the discovery of Higgs mass of around 
125 GeV in 2011-2012 and the consequent realization that a general framework such as 
the phenomenological MSSM (pMSSM) requires that the soft trilinear couplings A t ^ be 
large, we concluded [16] that the N MS GUT would be falsified by its failure to fit the down 
type quark masses unless Aq,[i were in the 10's of TeV : leading to a mostly decoupled 
superspectrum with only the LSP, gauginos and possibly a light slepton in the sub- TeV 
range ! The experimental data has now forced this realization on practitioners of MSSM 
parametrologv[T8] . In the NMSGUT it was a pre-diction. In the NMSGUT the successful 
fitting of fermion masses necessarily entails ultrasmall neutrino Majorana- Yukawa couplings 
leading to (first generation) right handed neutrino Majorana masses as small as 10 6 GeV. 
Taken together with the possibility of small values for the light generation Yukawa Dirac 
couplings it is possible to implement viable inflection point inflation by suitable tuning at 
the supersymmetric level itself. This is technically more appealing than a tuning applied 
to soft susy parameters which, being unprotected by SUSY, are unstable. We derive the 
tuning conditions for the NMSO(10)GUT and show how to satisfy them explicitly. 

In Section 2 we review and summarize the generic renormalizable single scalar inflaton 
inflection point model and calculate its slow-roll parameters, power spectrum and spectral 
index so as to use these results with the supersymmetric Type I Seesaw model(SIMSSM), 
with supergravity soft terms, once we have shown that it generates a suitable potential 
of the renormalizable type. In Section 3 we see how a generic Supersymmetric SU(3) x 
SU(2)l x U(1)r x U(1)b-l model, with generic supergravity type soft supersymmetry 
breaking terms, which contains the essentials of the Type I Susy seesaw implemented in the 
Minimal Susy LR Models and in R-parity preserving Susy GUTs, provides an attractive 
Inflationary scenario in parallel with its achievement of realistic neutrino masses. In Section 
4 we discuss the general features of reheating in this model and remark on the types of 
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Leptogenesis that can arise from the inflaton and right handed neutrino decay. In Section 
5 we give a discussion of the embedding in the NMSGUT and the inflationary parameters 
associated with realistic fits. We conclude with a brief discussion. 



2 Generic Renormalizable Inflection point inflation 

A generic renormalizable inflection point inflation model can be formulated in terms of a 
single complex field <p. Such a model[Sl[7] can reproduce the observed[3] inflationary power 
spectrum Pr = (2.43 ± 0.11) x 10~ 9 , spectral index n s = .967 ± 0.014 and scale invariance 
kdn s /dk ~ 0. It is also relevant to note that the ratio of tensor to scalar power spectrum 
amplitudes r = -p^ is known to be less than about 0.5. After extremizing with respect to 
the angular degree of freedom (which has positive curvature and cannot support inflation) 
one is left with the potential for a real degree of freedom (fr in the complex scalar inflaton 
field <p 



h 2 , 4 Ah ,o M 2 

— 4> 1=4> H 

VT 6^/3 2 



V = —<t ^ 3 + —<fr 2 (2.1) 



Here A, h, (fr are real and positive without loss of generality. The formulae we derive in this 
section are applicable to any single inflaton theory with a renormalizable potential. 

In the model of [B] the A, M receive dominant contributions from trilinear and quadratic 
soft supersymmetry breaking parameters: A, M ~ 10 2 — 10 4 GeV. A very small neutrino 
Yukawa coupling ~ 10~ 12 and a high degree of fine tuning between A and M is necessary to 
reproduce the observed inflation parameters [71 [8]. In our work however the contributions 
from soft supersymmetry breaking terms play a negligible role. The controlling mass scale 
is much higher , the required size of the yukawa couplings is larger and the degree of fine 
tuning is much less. 

It is convenient to trade the parameter A for a fine-tuning parameter A by replacing 
A = 4MV1 - A (A = j3 2 /4 in the notation of [7J). The inflection point at 

0o = ^(l-A + O(A 2 )) : V"(<f> ) = (2.2) 
is also a saddle point (V'((fro) = 0) when A = 0. For small A 

^o) = Vo = ^(l + 4A) ; V'(<t> ) = a = 
V"'(ct>o) = 7 = ^f(l-2A) (2.3) 

If the coupling h is tiny Vq » M 4 and 0o >> M/h. Notice that 7 tends to be quite 
small due to the smallness of h, while a is small (but non-zero[8]) because it is tuned to be 
small. The large vacuum energy and flatness of the potential around (fro then imply that 
if <fr starts with a value close to (fr = (fro and a small field velocity the universe will execute 
slow roll inflation as the field (fr rolls slowly down through a narrow field interval of width 
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A<p ~ Vo/jMp below <pQ. Around the inflection point <po, we can write the inflection point 
inflation potential in the form 

i 2 

V(<f>) = V + a{4> - <h) + \(.<t> ~ <k? + ~^{4> - 0o) 4 (2.4) 

The last term is essentially negligible since h 2 is very small by assumption. 
The slow roll parameters are defined as(M p = 2.43 x 10 18 GeV) 

M?V" M? 
ri{4>) = ~ -^7(0 - 00 ) 



£_ r-^ £_ 

^ 2 " ^ 



(2.5) 



The small first and third Taylor coefficients a, 7 determine [51 |2"0"1 [TU] the measured 
parameters of inflation (PR,n s ) once the field values {(^CMB^end) at the time of horizon 
entry of the "pivot" momentum scale (k p i vo t = 0.002 Mpc -1 ) and at termination of the slow 
roll are fixed [3 [19] (on the basis of an overall cosmogonic scenario and the consistency of the 
slow roll approximation (rj(4>end) ~ 1) respectively). k p i vo t corresponds to a representative 
scale of current cosmological observations. The field value at the beginning of inflation is of 
notional interest only. It is the number (Ncmb = N(4>cmb)) of e-folds of inflation left to 
occur after field value ^CMB reached at the time when the fluctuation scale of interest (k P i VO t) 
left the comoving horizon ( i.e k = akH^) during inflation that is of significance. This 
number is determined by the overall history of the Universe from primordial times [19|. 
Plausible inflationary cosmogonies require 40 < Ncmb < 60 and this severely restricts the 
inflation exponents. 

The field value at the end of slow roll inflation <p e nd is defined as the value where 

7 (</w ~ 0o)M 2 
V{<f>o) 

which gives 

V 



'end 



(2.7) 



Then in the slow roll approx <p ~ 0, <fi = —V'((j))/3H, where H = W V(4>o)/(3M 2 ) is the 
(constant) inflation rate during slow roll inflation, one has 



N(<f>) 



H 2 



V 

TV 



07 



m2 v arctan J — (</> - (f) end ) - arctan J — (<j) - <f))) (2. 
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and conversely 



<f>(N) 



<f>end + 00 (00 - 4>end)J ^ tan J 



1 + (00 - e „d) A /2^tan 




ajNMl 
2 V 



(2.9) 



It is worth remarking that this inversion of the function iV(0) was derived without 
assuming that (f>end 

« 4>(N) [5J. Together with an interpolating function derived below 
it allows us to obtain analytic formulae for the relations required among the parameters 
of the inflationary potential for successful inflation: avoiding tedious and opaque graphical 
methods [7]. 

The observed Cosmic Microwave Background(CMB) data [3] pose constraints on the 
power spectrum and spectral index for modes around the pivot scale. Barring non-standard 
scenarios where the post-inflationary period is punctuated by episodes of modified expan- 
sion, the number of remaining e-folds at the time the pivot scale left the horizon during 
inflation may be estimated by using the standard Big Bang thermal cosmogony along with 
estimates of the reheating behaviour of the universe after inflation. This gives[19j 



where p r h is the energy density after reheating and Vq the potential value during inflation. 
The reheating behaviour of the NLH flat direction inflaton is quite different from the the 
original quadratic chaotic inflation models for conjugate sneutrino inflation [21} I22j . The 
tripartite composition of the inflaton out of L, H, v c degrees of freedom will ensure that the 
bulk of the energy in the inflaton will be dumped into light degrees of freedom on the very 
first oscillation. For the present we merely assume that the reheating is immediate so that 
one can set p r h = Vq in Npi VO f. We then find that for M in the range 10 6 — 10 12 GeV, h lies 
between 10 -95 to 10 -65 and then N p i vot = Ncmb = 51 ± 5 adequately covers the possible 
range. Even if this range is lowered by effects of reheating or non-standard cosmogonies the 
effect on the relevant exponents will prove to be marginal. Although the observed CMB 
is actually a combined spectrum of modes exiting the horizon around \N — iVpi vo t| < 5, 
we can approximate and regard it as the single spectrum from the mode that exits the co- 
moving horizon when = (ftcMB only. Thus = 4>cmb is the field value near 0o where the 
inflation giving rise to observable effects today kicks in (when Ncmb e-folds of inflation are 
remaining). The power spectrum and spectral index we see today are then -Pr(0(-/V"cmb)) 
and n s (0(iVcMB)) respectively, where \N C mb - N p - lwot \ < 5. 

The slow roll inflation formula for the power spectrum of the mode that is leaving 
horizon when the inflaton rolls to is ([20 ) 



= 65.5 + In 




(2.10) 



pivot 



M p 



p r {4>) 



V 



(2.11) 



247r 2 M4 e (0)' 
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and the corresponding spectral index and it's variation with momentum is 



n s {4>) = 1 + 2ry(0) - 6e(<£) 
Vk (n s ) = ^M = _l 6e?? + 246 2 + 2e 2 (2.12) 

The ratio of tensor to scalar perturbations r = -pj^ = 16e. In practice e,£ are so small in 
the narrow region near (f)Q where slow-roll inflation occurs that their contribution to n s is 
negligible. Thus T>k(n s ) is negligible i.e. the spectral index is scale invariant in the observed 
range, as is allowed by observation so far. 

To search for sets of potential parameters M,h,A compatible with Pr,us, Nqmb in 
their allowed ranges one may proceed as follows. First one uses the chosen (within experi- 
mental range) values of Pr, n s and given M, h to define 

V 



(-CMB 



24m 2 M$P R 



VCMB = 5 ( 2 - 13 ) 



From e c mb, Vcmb one may deduce o-cmbAcmb using the eqns.([23 

Vq^cmb 



4>CMB = 4>0 + 



jM 2 



V 

2^2 



acMS = (2.14) 

The required fine-tuning A is then 

A = = (-^r)\^^= " (1 " n s f) (2.15) 

otcMBi^- should emerge real and positive and using {ocmb, 4>CMb} in the formula for 
Ncmb one should obtain a sensible value in the range Nqmb = 51 ± 5. Positivity of A (a 
local minimum develops if A is negative leading to eternal inflation) requires 

9 3vr(l -n s ) 2 JP~R~ 

h 2 > M — v (2.16) 

16M P v ; 



Using eqns. (|2.3|2.7tt2. 13j) in eqn (|2.8|) we have 

1 2z(l + n 

— arctan — ^ 

z 8z z + 1 — n 



Ncmb = - arctan - _^ ■ - — ^— (2-17) 



( h2m > - (i ~:- )2 )i (2.i8) 



^kM^/PIi 16 

By solving eqn. (|2.17p for z = zo(NcMB,n s ) one obtains the general relation between h and 
M : 

h 2 3ttJP~r, 9 ,„ t , (l-n s ) 2 N 

M = -^^(N C MB,n s)+ '— T ^-) (2.19) 
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and then 



16M 2 z 2 (N C mb, 



97T 2 M 2 P R ((l-n s ) 2 + 16z 2 (N C MB,n s )) 2 K ' ' 

Where zo(NcMB,n s ) is the solution of eqn (|2.17|) . An excellent approximation to the the 
required function in the region of interest in the Nqm b > n s plane is given by the Taylor 
scries around n° s = 0.967, iVg = 50.006 : 

z (Ncmb, n s ) = .0238 - 0.0006(iVc7MB " N%) + 0.2407(n s - n° s ) 

{n C mb-n° c ) 2 _ q 7nS 7^K-"g) 2 xnnn^Af^,,. _ /vQ.V n _ n o^ 



+0.000022 liVCMB 2 1 " c> - 3.70875 ( " s ~ nJ + 0.002353(A^ M B - N°)( 



_ m00015 (JfcMB_N^ + 8 79982 (n £ ^)! _ 0,000788 (AfcMa - jV f )2(ns " n ° ) 

-0.7536 (jVcMB ' JV | )(Wa " w ° )2 (2.21) 

In Fig. 1 we have plotted the contours of zo{NcMB,n s ) in the NcMB,n- s plane and one 
sees that the variation of zq is rather modest. So for the plausible range 46 < Nqnb < 56 
one obtains a tight constraint on the exponents in the relation between h,A and M: 

1,2 -in~24.95ifc0.17,' \ . a i n~28.17±.13 ( j*j \2 fo oo\ 

h ~ 10 (_) , A ~ 10 (_) (2.22) 



We have estimated the maximum variations in the exponents corresponding to the 
quoted errors in the WMAP 7- year data [3 J from the graphs in Fig. 2 and Fig. 3. 

However a clearer qualitative understanding results from noticing that for Ncmb ~ 50, 
Zq rs N ^ IB solves eqn. ([2.17p to a good approximation. Then eqn. (|2,18p gives 

h 2 _ 3vr v^Pr 2.75 x 10~ 22 _ in _ 25 



Af ~ Mp N^, ArR ~ iV 2 M ~ 10 ~ ( 2 - 23 ) 

10" 28 - 2 Gey^ 2 (2.24) 



"CMB CMB 

A 4.14 x 10~ 34 



^ 2 ^ 

Thus these simple approximate expressions give effectively the same results as the more 
carefully derived expressions in eqn. ()2.22p . Thus we have viable inflation with 

M 4 

V ~ — ~ (M) 3 x 10 25 GeV ~ 10 43 - 10 61 GeF 4 (2.25) 
n/ 



1 



J^~10 3 -10 12 Gey ; T max ~ V 3 ~ 10 11 - 10 15 GeV (2.26) 



It is clear from eqn, (|2.22p that the fine-tuning measure grows with M so that /3 = vA 
can be as large as 10~ 2 for M ~ 10 12 GeV. Due to the large value of the inflaton mass 
compared to the case of MSSM inflation[4j or Dirac neutrino inflation^ [8] the fine-tuning 
of parameters required is much less severe and no additional dynamics need be invoked 



9 




0.955 0.960 0.965 0.970 0.975 0.980 



Figure 1: zq contours in the NcMB,n s plane. The variation shown contributes to the small 
range of permitted magnitudes for h 2 /M,A/M 2 etc 
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Figure 2: Variation of exponent of jj with Nqmb for different values of ti s ,Pr. 
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Figure 3: Variation of exponent of ™ with Ncmb for different values of n s; Pr. 
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to make it plausible[71 [23]. It is also important to note that the ratio of tensor to scalar 
perturbations r ~ 16e ~ 2(M/10 14 GeV) 3 . Since M is at most the heaviest right handed 
neutrino mass ~ 10 12 GeV it is clear that it is difficult to get r > 10 -6 . Thus measurement 
of tensor perturbations via the Cosmic Microwave background polarization at the r ~ 10~ 3 
level or larger would not be compatible with inflection point inflation controlled by the right 
handed neutrino mass. Any renormalizable single inflaton model must respect these generic 
constraints and yield values of its associated parameters that are sensible in terms of the 
other (particle) physics that it describes. It remains to specify the Renormalizable Susy 
seesaw Inflaton scenario and consider the NMSGUT as a self contained realistic test bed. 

3 Supersymmetric seesaw Inflaton model 

The essentials of the Supersymmetric seesaw inflation scenario may be captured by consid- 
ering a model with gauge group SU(S) x SU(2) x U(1)r x U{1)b-l and the field content 
of the MSSM with some additional superfields. Soft supersymmetry breaking terms are of 
the supergravity type [i.e trilinears proportional to yukawa couplings and universal, or uni- 
versal except for Higgs (NUHM scenario), soft scalar masses]. The essential fields beyond 
the MSSM consist of a right handed Neutrino chiral multiplet N[l, 1, —1/2, 1] and a field 
S[l, 1, 1, —2] whose vev generates the large Majorana masses M v (10 6 — 10 14 GeV) for the 
conjugate neutrinos v c A = Na via a renormalizable superpotential coupling 3\/2/ ' abSv c a v c b . 
Additional fields Gj which serve to fix the vev of S are also present as in Minimal Super- 
symmetric Left Right Models (MSLRMs)[l3] and in GUTs that embed them [Hl fTo l fTB]. 
The other essential component of the scenario is neutrino Dirac mass generating Yukawa 
couplings yAB, A, B = 1, 2, 3 in the superpotential. These couple the right handed neutrinos 
to the Left chiral lepton doublets La = [y e) A , A = 1, 2, 3. La transform as L[l, 2, 0, —1] 
and the up type Higgs doublet type field as H[l,2, 1/2,0] so that tjabNLaHb is a gauge 
invariant term in the Superpotential. Of course each such doublet present in the underlying 
theory must have its complementary doublet transforming as e.g. H[l, 2,-1/2,0] to cancel 
anomalies. The relevant flat direction is assumed to extend out of the minimum of the 
supersymmetric potential corresponding to the breaking of the gauge group down to the 
MSSM symmetry 

SU(3) x SU(2) x U(l) R x U{\) B -l -> SU{3) x SU{2) x U(1) Y (3.1) 

This leads to a Type I seesaw plus MSSM (SIMSSM) effective theory. 

After the breaking one has Y = 2T^b + {B — L) where T^r is the U(1)r generator. Note 
that unlike the case of the Dirac neutrino masses scenario [B] B — L is not a gauge symmetry 
down to low energies. This can have important consequences for nucleosynthesis and matter 
domination since the heavy right handed neutrinos must find a non-gauge channel to decay 
through. In the present case this channel must perforce be a Yukawa coupling since the 
right-handed neutrinos are singlets of the low energy (SM) gauge group. This is in contrast 
to the Dirac scenario where a low scale of B-L breaking is assumed so that v c can decay via 
gauge couplings. 
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The fields TV (the chosen conjugate sneutrino ), v{ chosen left sneutrino flavour from a 
Lepton doublet L with suitable Yukawa couplings) and the light neutral Higgs field (from 
the doublet H with Y = +1) may be parametrized in terms of the flat-direction associated 
with the gauge invariant NLH as 

N = t> = h = ^= = <fte ie ; (ft > 0, € [0, 2vr) (3.2) 
v3 

The additional fields Oj, unspecified at the moment, are assumed to be coupled to S in such 
a way that extremization of the SUSY potential using = 0, -D a |0=o = fixes the vev 
of S: < S >= a/y/2 without constraining the inflaton field (p. This is of course true in the 
Minimal Susy LR models[13j and renormalizable Susy SO(10) GUTs [151 HE] which are our 
inspiration. 

The vanishing of the D-teim for the B—L generator requires Vti to include the companion 
field(s) S[l, 1, —1, 2] which have a vev of equal magnitude as S in order to preserve SUSY 
through the symmetry breaking down to the MSSM symmetry at high scales. This is just 
as in MSLRMs and R-parity preserving GUTs |13t [141 [i~5| 116] . The gauge invariance of 
NLH ensures that the D-terms for the flat direction vanish. Thus at scales (ft ~ a >> Mg 
where SUSY is exact the relevant superpotential is given by: 

W = 3V3yNuh + 3fV2SNN + ... = yp 3 + fV2Stp 2 + ... (3.3) 

where h, f, a can be taken real without loss of generality. The right handed neutrino Ma- 
jorana mass will be M v = 6 fa. 

Since the equations of motion of the unperturbed vacuum imply < F$ >= 0, 
< S >= a/y/2 this superpotential leads to a flat direction potential 

V SUS y = |3^ 2 + 2M9| 2 + 2|jV 2 | 2 

= / 2 [(2 + 9y 2 )0 4 + I2y<p 3 a cos 9 + 4ct V] (3.4) 

Here y = y/f and we see that a sets the mass scale. Minimizing with respect to 9 gives 
9 = ir. In so far as we are here interested only in the inflationary dynamics (once parameters 
have been tuned to ensure an inflection point in the plateau region where \ip\ ~ a) we can 
focus on just the real part of (p and set tp = —(ft with (ft real and positive near the inflection 
point but free to fall into the well around (ft = and oscillate around that value. The 
imaginary part (ft' of ip has a large curvature Vgy ~ <r 2 in the plateau region. Since 
V0'|^'=o = it is consistent to consider the dynamics in the real p plane alone as a leading 
approximation. The effect of jitter in the (ft 1 direction when the dynamics is initiated with 
(ft' can be studied numerically as a correction to the dynamics of the inflaton field (ft. 

In addition one also expects a contribution to the potential from the /i term for the 
Higgs doublets together with SUSY breaking quadratic and cubic soft terms, which we 
assume to be of the type generated by supergravity, but with non universal Higgs masses, 
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i.e of the form: 



Vsoft = [Mw 3 + /V2V) + h.c] + m 2 f J2\f\ 2 + m n\ H \ 2 + m W\ 2 

f 

= f 2 \ y A (/) 3 a cos 39 + A a 2 <p 2 cos 29 + mla 2 (I) 2 ] (3.5) 

here rho = mo/ fa, Aq = 2Aq/ fa. The soft mass mo receives contributions from the sfermion 
and Higgs soft masses as well as the /x term 

2m 2 - + m 2 H 

m 2 = (3.6) 

mj H are the sfermion and up type Higgs soft effective masses at the unification scale 
(jn 2 H = m 2 H + |/i| 2 ). Since these masses and Aq should be in the range 10 2 — 10 5 GeV while 
the righthanded neutrino masses lie in the range 10 6 — 10 12 GeV, it is clear that rho, Ao are 
small parameters and even for the large values of mo, Ao ~ 10 5 GeV found in the NMSGUT 
fho,Ao « 1. Thus these terms cannot significantly change 6 = it assumed earlier. The 
total inflaton potential is then 

Vtot = f ((2 + 9y V - (Ao + 12)ya^ 3 + (i + m 2 + A)a 2 <p 2 ) . (3.7) 

Thus we have a generic quartic inflaton potential of the same type as in Section 2 but 
the parameter values in the case of Type I seesaw are quite different from the light Dirac 
neutrino case. We have the identification of parameters 



h = iVl2(2 + 9y 2 ) 
A 3/(i + 12)ya 

M 2 = 2f 2 a 2 (4 + A + m 2 ) 
A 2 , 



A 



16M 2J 

9y 2 (A + 12) 2 
32(2 + 9y 2 )(I + " : i 2 +4) 



(3.8) 



For seesaw models the natural magnitude for the neutrino Dirac mass is, mj/ > lMeV 
(i.e \y„\ > and then the limit m v << O.Oley for the lightest neutrino (assuming direct 
hierarchy) implies M v c > 10 6 GeV). Since the preferred values for the Susy breaking scale 
are smaller than 100 TeV (at most) it follows that the maximum value of |-Ao|> \ mo\ ~ 0.1 
and they could be much smaller for more typical larger values of the conjugate neutrino 
masses M„c ~ 10 8 to 10 12 GeV. It is then clear from the corresponding range A ~ 10~ 12 to 
10 -4 that the coupling ratio y = y/f becomes ever closer to exactly y = 4/3 as M increases 
and even for M ~ 10 6 GeV differs from 1.333 only at the second decimal place. Thus to a 
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good approximation h = 6y6f. Then it follows from the Eqs. (|2.22p and ()3.8[) that 



f „ ln -26.83±0.17/ a \ . Ti/r ^ -ir.-25.38zfc0.17/ a n2 

J GeV ' GeV 

A ~ 10 -78.93±0.47 ( ^ ) 4 ,3 Q) 

GeV 

The range M ~ 10 66 to 10 10 6 GeV corresponds nicely to 10 16 GeV < a < 10 18 GeF: as 
is natural in single scale Susy SO(10) GUTsJTH [151 EH E|- / increases with a with values 
below to 1CT 11 achievable in the NMSGUT only with difficulty. Of course in MSLRMs, 
since there are no GUT constraints on a, one can assume somewhat wider ranges for these 
parameters. 

In all relevant cases A < 10 -4 is required. Thus the above equations imply that y 2 must 
be close to the value 

9 16 - 8^0 - 32m 2 + A 

Here Ao,mo ~ 0(M$ /M u c) << 1, hence j/o is rather close to 4/3 and the equality is very 
close for larger M ~ fa since then Aq,itio are tiny. This then is the type of fine tuning 
that supports the development of inflation in SIMSSM models. We see that the measure of 
severity of fine tuning (3 = \/A ~ 10 -2 — 10~ 6 compares quite favourably with the case of 
the MSSM or Dirac neutrino inflaton since there (3 ~ 1CP 12 to 10~ 10 due to the low values 
of the inflaton mass in those cases. The dominant component of the fine tuning in the 
present case is a fine-tuning of superpotential parameters, which is radiatively stable due 
to non renormalization theorems. Specially for large a > 10 16 GeV the Type I Susy seesaw 
can provide a rather attractive inflationary seesaw with a natural explanation for neutrino 
masses and weaker tuning demands on the radiatively unstable Susy breaking parameters 
than the extreme and unstable fine-tunings demanded by typical inflection point scenarios 
and in particular the Dirac neutrino model [6]. Moreover, unlike the chaotic sneutrino 
inflaton scenario [21] [22], no trans-Planckian vevs are invoked. 



4 Reheating and Leptogenesis 

After inflation concludes the energy stored in the inflaton will be transferred into a thermal 
bath of the MSSM degrees of freedom. Determination of the time required to thermalize 
the inflaton energy and the resulting reheat temperature T r /j (i.e the maximum tempera- 
ture of the thermal bath after thermalization) requires understanding the post-inflationary 
dynamics of the LHN flat direction inflaton. An important issue that can be tackled at 
the level of the effective SIMSSM is generation of the the cosmological baryon number 
asymmetry(ns/n 7 ) via Leptogenesis [lOj . Although a detailed analysis of these issues re- 
quires a separate publication, the existence of previous detailed studies of preheating [TT] in 
a MSSM flat direction inflaton model[24j and of non-thermal Leptogenesis in a preheating 
model[12j make the generalizations required to combine the two ideas in the context of su- 
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persymmetric seesaw inflation easy to outline, but too long to derive, here. Super symmetric 
seesaw inflation offers an attractive synthesis fulfilling the need expressed in [12! : 

" There have been many models of leptogenesis. A hallmark of our model is the economy 
of fields. The only undiscovered fields are the inflaton, (f>, the standard model Higgs, h, and 
the right-handed neutrino, N. There are very good reasons for suspecting that all exist! The 
only unfamiliar aspect of our model is the strong coupling of the inflaton field to the Higgs 
field. While there is no reason to preclude such a coupling, it would be very interesting to 
find particle-physics models with a motivation for the coupling. " 

In our model the the inflaton is itself partly comprised of the Higgs field and therefore 
fulfills the requirements of |12j exactly, besides bringing together a number of other related 
streams of thought. We remark however that the situation is made more complex by the 
high reheating temperature associated with the large inflaton mass. Thus both thermal and 
non-thermal leptogenesis may contribute to the generation of n^/n 7 . 

Due to the gauge(H,L) and third generation yukawa(-ff) coupled components of the in- 
flaton the inflaton energy is likely to decay very rapidly (i.e in decay time « ~ 
(hM p ) I M 2 ) through the so called 'instant preheating" mechanism [1 1 \ fl"2| I24j . In this mech- 
anism the preheating dynamics results in a rapid decay (well within a Hubble time) of the 
complete inflaton vacuum energy into a radiation bath which therefore thermalizes to a 
temperature determined essentially by the equality between the radiation bath energy and 
the starting inflaton energy. This gives an estimate for the reheating temperature 

T rh ~ T max ~ y 1/4 ~ M/h 1 ' 2 ~ 10 11 - 10 15 GeV (4.1) 

The parametric dependence is identical to that found in [23], the difference in scales arises 
only because the inflaton mass M ~ 10 6 — 10 12 GeV in our model is much larger than the 
inflaton mass ~ 0.1 — 10 TeV in [21] coming from soft Supersymmetry breaking. 

In the preheating mechanism a class ( "x type" ) of degrees of freedom, whose masses(m x ~ 
g4>{t)) and decay rates (T ~ g s cp(t)) are proportional to the instantaneous inflaton value 
4>(t), are produced non-perturbatively every time the instanton field crosses zero. This oc- 
curs since the x modes are ultralight for a sufficiently large time interval around the zero 
crossing time during which adiabaticity is violated ( Co}, > oj 2 : where is the oscilla- 
tion frequency at wave number k). In our model the x modes are the components of the 
H, L, u1,ul chiral superfields and the W±, B gauge superfields. In fact the x modes can be 
identified simply by checking which fields become massive in the presence of background val- 
ues of the three components of the inflaton [y, z/£, h°). Then with the usual superpotential 
(we have suppressed generation indices) 

W = y u Q L Hu c L + y d Q L Hd c L + y u LHN + y l LHe c L + ... (4.2) 

we see that y u leads to massive ul,u^, y u leads to massive (one combination of the 
three cl ), h°, h + , vl, y l leads to massive h~~,e c L (one combination). Since < H, N, L > 
preserve U(l) em , the gauge couplings give masses to Z (which forms a Dirac supermultiplet 
with {y — ho)/y/2) and W± (form a pair of Dirac supermultiplets with l~,h+). This set 
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of fields are the x tyP e fields whose mass varies strongly with as it oscillates and whose 
production, when ~ 0, and decay when eft >> M\y, is the basis of 'instant preheating'. The 
inflaton vev leaves the down quark and gluon/gluino fields and ho, and some combinations 
of the 17,1% fields with light (MSSM type) masses. These light (f/'-type) fields will form 
the first step in the decays of the x field. As < <fi > again increases the x modes become 
very heavy and unstable and as a result decay rapidly(within a time T& ec ~ « m^ 1 ) 
to the light (mostly coloured) MSSM d.o.f. to which they are coupled (dominantly via the 
D-terms and gauge-yukawa terms but also via the Superpotential couplings for the third 
generation). As a result a significant fraction ~ 1CP 1 of the inflaton condensate energy 
passes into the light MSSM modes with every crossing resulting in complete transfer within 
~ 10 2 oscillation times. 

rose ~ m^ 1 « flT* Jn ~ (hM p )T osc /M ~ (1 - 150)r osc (4.3) 

Once the energy is in the light modes MSSM interactions, in particular the gauge interac- 
tions, are sufficient to rapidly complete thermalization so that essentially all the inflaton 
energy will be thermalized within, at most, a few Hubble times after the end of inflation. 
Rapid decay of the inflaton oscillation amplitude leaves the light modes to thermalize the 
energy dumped by the inflaton into a radiation bath of all modes: which are no longer ever 
heavy because the inflaton has decayed. The reheating temperature is 

T rh ~ (4^) V V 4 ~ T max ~ 10 11 - 10 15 GeV (4.4) 

where g* = 228.75 is the effective number of MSSM degrees of freedom. The essential 
point is that the reheating temperature is well above that required to produce relativistic 
populations of gravitinos : which are unacceptable if their lifetimes are larger than the nucle- 
osynthesis time tat ~ 1 sec since their decay after nucleosynthesis would destroy the created 
nucleons. The straightforward and generic resolution of this gravitino problem is if the 
graviton masses are sufficiently large so that the gravitinos decay before nucleosynthesis [25] 

Tqrav ~ 10 5 Sec( ) 3 << T N ~ lseC (4.5) 

Thus we see that the viability of Super symmetric seesaw Inflation strongly indicates that 
the scale of supersymmetry breaking -as indicated by the gravitino mass- should be above 
50 TeV. The fact [16] that such large supersymmetry breaking scales are preferred by both 
the NMSGUT and the latest data indicating[26] light Higgs mass Mh ~ 125 GeV rounds off 
the picture nicely. Furthermore such large reheat temperatures also ensure abundant ther- 
mal production of all flavours of righthanded neutrinos after inflation. Their CP violating 
decays into leptons can generate the net lepton number density which drives creation of the 
requires n#/n 7 by Sphaleron processingjTD] . Thus the NMSGUT can not only accommo- 
date inflation but is also compatible [27] with (thermal) Leptogenesis [10] for generating the 
observed baryon to entropy density ns/s ~ 10~ 10 . 
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An interesting additional source of non-thermal leptogenesis is provided when one 
realizes [T2] that, the Higgs field H, which is one of three fields making up the putative 
inflaton direction in field space, is itself a x tyP e field and is furthermore coupled to the 
righthanded neutrinos. Thus during the course of oscillations of the field components of the 
inflaton - which commence at the end of the inflation - the Higgs mass mh ~ g2<ft fluctu- 
ates to values both below and above the righthanded neutrino masses which are essentially 
constant at M u c ~ fa even in the presence of the inflaton (i.e D, h, N) background since 
§2 » f,y- Thus as the Higgs mass oscillates below and above the right handed neutrino 
masses one expects CP violating -therefore net lepton number producing - inter-conversion 
of the Higgs with righthanded Neutrinos as in [12] . If this net lepton number is not washed 
out by the inflaton energy dump (so that a Hubble volume contains a certain Lepton excess 
produced by this inter-conversion even though average energies are well above the mass of 
the right handed Higgs) then we may expect that a non-thermal Leptogenesis component 
will add to the thermal leptogenesis due to decay of the righthanded neutrino bath. 

An important complication in the present case, that we have glossed over in the above 
account, is that the L, H and N components of the inflaton can have quite different decay 
rates once the gauge interactions are effective, since iV is a gauge singlet. A proper analysis 
must track the evolution of all three fields making up the inflaton - from an initial condi- 
tional^ end of inflation) where they start out equal. This makes the equation of motion 
and Boltzmann equation for the relevant degrees of freedom significantly more complex 
and this requires a separate numerical study which involves the interplay of the couplings 
§ A i Dab 1 92 ■ The study of this evolution and the operation of Leptogenesis in these models 
is now in progress. 

5 Inflation and neutrino masses in the NMSGUT 

Finally we consider the embedding of our generic Type I scenario in a realistic Susy SO(10) 
model[16"l [T71 [28] that has successfully fitted the known fermion mass-mixing data and can 
also be consistent with limits from B violation and other exotic processes [17j. We will see 
that neutrino flavour plays a key role in enabling inflation : the model favours an inflaton 
composed of third generation conjugate sneutrino, first generation left slepton (sneutrino) 

and T m = 1/2 Higgs. 

The New Minimal SO(10) GUT (NMSO(IO)GUT) uses Higgs fields in the 210, 126, 126 
representations of SO(10) which contain 5 SM singlets whose vevs break SO(10) down to the 
SM gauge group at a superheavy scale Mx- Three of these vevs, called p, u, a come from the 
210-plet and one each from the 126(<r), 126(a). An explicit solution to symmetry breaking 
50(10) — > SIMSSM, in terms of a simple cubic equation for a complex variable x and 
depending on a single parameter ratio £; was found in the third paper in [15]. This solution 
preserves supersymmetry and makes no use of the soft breaking terms which constitute a 
negligible perturbation of the global susy symmetry breaking problem [29], in the sense that 
they modify the superheavy vevs ~ M x ~ 10 17 GeV only by terms of order M s ~ W 4 GeV. 
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The spectra calculated |1 5 1 [T6] using this analytic solution for the the MSGUT vacuum 
are the basis of our detailed Renormalization Group analysis of grand unification in this 
class of models. Questioning the received wisdom that large SO(10) representations make 
grand unification futile [3"U] we showed [TBI El EI] that the inclusion of threshold corrections 
considerably ameliorates the problem of large gauge beta functions by allowing one to raise 
the threshold corrected unification scale close to the Planck scale and lower the gauge 
coupling at unification. Taken together these features imply that even with the huge beta 
functions characteristic of MSGUTs the problem of a Landau pole in the gauge coupling 
may be postponed to the Planck scale : where it becomes moot along with the structure 
of space time anyway. The physics of asymptotically strongly coupled gravity and gauge 
theories is anybody's guess (see however [32J for our speculations and simplified model for 
'tamed' asymptotically strong GUTs). There are even claims that gravity is capable of 
ensuring the asymptotic freedom of any gauge theory [331 ;Mj. It is also possible that a RG 
fixed surface on which the gauge coupling remains weak in the UV may exist. In view of 
the many uncertainties we take the stand that the large beta functions of the NMSGUT 
are not an issue that need prohibit the study of these minimal and realistic theories. 

The grand unified minimum of the potential defined by the vevs f2 = {p, u, a, a, a} 
shifts only by fractions of order 10 -12 due to supergravity mediated soft supersymmetry 
breaking terms. The D terms of SO(10) are all exactly zero for these vevs. To examine 
the issue of an inflaton corresponding to the NLH flat direction in the SIMSSM we must 
demonstrate the existence of a corresponding flat direction of the full GUT potential based 
on light (SIMSSM) field vevs. This flat direction rolls out of the grand unified minimum that 
defines the MSGUT vacuum with the SIMSSM as its effective theory. The relevant fields 
are the GUT scale vev fields f2 and the (6) possible components hi, h{;i = 1...6 of the light 
MSSM Higgs doublet pair H, H together with the chiral lepton fields L A ,v c A ,A = 1,2,3. 
The relevant superpotential is then[T5J [TE] 

W = 2\f2{h AB h x - 2V3f AB h 2 - g AB (h 5 + iV3h 6 )) + h T H(< Q >)h 

+AV2f AB d9 A v B + W n (0) (5.1) 

where 

Wa{n) = m(p 2 + 3a 2 + 6w 2 ) + 2A(a 3 + 3poo 2 ) 

+ (M + r](p + 3a - 6u))aa (5.2) 

and 

d Wo 

-Q^\h,u,L=o = D a (n)\ h>P , L=0 = Q (5.3) 

here h AB , g AB , f AB are the yukawa coupling matrices of the three matter 16-plets to the 



10, 120, 126 Higgs multiplets respectively. Equation (|5.3p defines the MSGUT vacuum[15J. 

Of the 5 diagonal D-terms of SO(10) only those corresponding to the generators T^l, 
T3R, B — L are charge and color neutral with vevs for Q and v, u c , h^. The vevs ft do not 
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contribute to these D terms so their values are 

^3L = f (-X>o| 2 + £M 2 ) (5-4) 

i=l A 

D *R = y (E l^ol 2 " 2 l^ol 2 - £ l^l 2 ) ( 5 - 5 ) 

Db—l = y|5«(E(l^l - l^l 2 ) + 2 l /l4 °| 2 ) ( 5 - 6 ) 

where we have used the fact (paper 5 in [15] and [16]) that only h± a = $| has B — L = 
-\-2,T^pi = —1/2 and thus T^l = 1/2 while all others have T^r = 1/2 and B — L = 0. 
Note that g u is the SO(10) gauge coupling in the standard unitary normalization. Thus the 
D-flatness conditions are 

m 2 = E i^oi = E M + 2 im 2 ( 5 - 7 ) 

For simplicity we assume that only one generation each of sneutrinos va and conjugate 
sneutrinos 9b contributes to the inflaton flat direction; but not that they must belong 
to the same generation. At this point we remind the reader [T5| [T6] that in MSGUTs 
the MSSM Higgs doublet pair is defined by fine tuning Det{%) ~ so that its lightest 
eigenvalue [i ~ Mjy ~ 1 TeV specifies the [i term in the superpotential of the SIMSSM : 
W = fiHH + .... The doublet pair H, H is a linear combination [T5"l Ell ES] of the 6 doublet 
pairs of the the NMSGUT : 

hi = UijHj hi = UijHj (5.8) 

where U, U are the unitary matrices that diagonalize the doublet mass matrix % : UT-LU = 
Diag{fi, M% , Mq 1 } to positive masses. To leading approximation they can be calculated 
with [i = = Det(H). The so called Higgs fractions : on = Un,cti = Un , are crucial 
in determining the grand unified formulae |15[ [T6] for the fermion yukawa couplings that 
give rise to the fermion masses. To obtain the tree level yukawa couplings one makes the 
replacement hi, hi —¥ ctiH,aiH in the expressions coupling the GUT Higgs doublets(/tj, hi) 
to the matter fermions of the SIMSSM. Thus in particular the neutrino Dirac coupling is 
({hAB,gAB,fAB)=^V^(hAB,9AB,fAB) ) 

Vab = hABOL\ - 2V3/abU2 - 9As(a5 + iV3at e ) (5.9) 

From the |-FX| 2 contributions to the potential it is clear that the involvement of any but the 
light Higgs doublet H would lead to GUT scale rather than conjugate neutrino scale masses 
for the inflaton. Moreover in view of the stringent upper bounds on the fermion yukawas 
(see eqn. (|2.22p ) the involvement of the lightest generation is unavoidable. Thus we take 
i>A = v\. However if we also take va = v\ we find that the tuning constraint has at best the 
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form 1 2/1 1 1 2 ~ 10(| J/21 1 2 + | J/31 1 2 ) : which is very hard to satisfy with normal neutrino mass 
hierarchy(the case studied so far) in (N)MSGUTs . On the other hand with v A = v% there 
is a possibility of satisfying the fine tuning condition. Thus our ansatz for the flat direction 
fields is 

z^i = — h i0 = — v 3 = — VI - 2\a 4 \ 2 (5.10) 

Notice the peculiar role of the Higgs fraction 04 which enters the flat direction ansatz as as 
r = 1 — 2|<3!4| 2 . As it happens the solutions we have found earlier [16] often have {a^l ~ 0.5. 
Thus it is not inconceivable that T can be consistently tuned to zero by varying the GUT 
parameters. The challenge is to do so without destroying the realistic fits to the fermion 
data. 

By varying the fields Q,,va,va we can now easily derive the F-term potential 

v h ard = [(yV)n + r(|/i3i| 2 + 4|93i| 2 + (yV t )33) + 4|/33l 2 r 2 ]^ 

+^h3\y u 3 1 \\v\VTcos(e^ + 6 y » i -0*)M 3 + (5.11) 

[^ + fl/ 3 3lVl 2 r]M 2 (5.12) 

We can also write down the generic Supergravity(SUGRY)-NUHM generated soft terms in 
terms of a common trilinear parameter Aq but different soft mass parameters m 2 ~,mf t _ for 
the 16 plets and the different Higgs (we have dropped the constant term from M$W{Q) 
assuming it is removed by the Supergravity scenario tuning to set the GUT scale vacuum 
energy to zero by tuning hidden sector parameters). The differences among the SO(10) 
Higgs soft masses could be due to renormalization from the threshold corrected unification 
scale/Planck scale to the scale M° x = 10 1625 GeV at which the SIMSSM and NMSGUT are 
matched in our work |36j. 

Vsoft = A W + c.c. +m 2 6 |§| 2 + ^m|j/ii| 2 

i 

4 



2^0^1^1^=005(3^ + o y »j + -A,/ 3 3HrM 2 c 0S (2^ + e 9 ) 

+(a 2 -^)|0| 2 (5.13) 



where 



-9 mh, . ^fh? h \ai\ 2 \n\ 2 

™o = -yHi + r) + - Jh ^~ L + -Y- ( 5 - 14 ) 

i 

and rhiQjThhi, Aq are all ~ O(Ms) Now the extreme dominance ^33 1 cr | >> Ms implies that 
the phase 9$ is fixed by minimizing just the term in Vhard'- 

6<p = 7T + Q- a - Q y u (5.15) 
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We shall assume that 6$ is fixed at this value. Since the inflationary dynamics is at large 
values of \(j)\ and fixed 9^ we can work just with a real field 4>. Comparing the sum of the 
hard and soft potentials with the generic renormalizable inflaton potential in Section 2, we 
immediately obtain the parameter identifications 

h = ^[(yV)n + r(|/ l 3i| 2 + 4|53i| 2 + (yV t )33) + 4|/33| 2 r 2 )]^ 

A = i(16|/ 33 ||y|fi||ff|v^ + 4|^ 1 |v^^oCos(30 ff -2^ 1 )) 
M 2 = — \h 3 \ 2 \a\ 2 T + -A f 33 \a\TCos(39 & -28 y vJ + 2m 2 (5.16) 



The fine tuning condition A = 4M now becomes 



1 2/31 1 = 9A d -8A + '^ 21 ' +r d /l31 ' + 4 l33il + I2/32I +I2/33I ) 



+4i/ 33 rn ( 5 - 17 ) 



Where 



A ^ = ( 1 + ^ Cos (.^- 2B v'J) 2 ( 5 - 18 ) 



and M3 = /33|<t| Note that in view of the ratio between the soft breaking scale and the 
mass of the heaviest right handed neutrino, A nj d are both very close to unity. Thus the fine 
tuning condition is essentially between hard parameters as in GUTs and in sharp contrast 
to MSSM inflaton models @]: 

l2/3lP = I ^ f (r(|/ i 3l| 2 +4|53l| 2 + |^2| 2 + l^3| 2 ) 

+l2/nl 2 + l2/2il 2 + 4|/ 33 | 2 r 2 ) (5.19) 

In NMSGUT fits of the fermion data we typically find a strong hierarchy 1 2/33 1 >> | j/32 1 >> 
1 2/31 1 >> 1 2/21 1 > \yii\- So it is evident that one must tune 

T i.e |a 4 | « -7= (5.20) 

v2 

to a good accuracy. This means that the MSSM doublet H is almost exactly 50% derived 
from the doublet in the 210 plet ! If this condition can be achieved the remaining tuning 
condition is only 

I2/31I 2 = 8(|tf 1 | 2 + |tf 1 | a ) (5.21) 

which is easy to enforce in the NMSGUT. 
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However there is an additional demand coming from eqn (|2,22p : h 2 /M3 ~ (y v 'y v )ii/Ms ~ 
10 -25 which is, at first glance, much harder to enforce. It is rather remarkable that our 
results in[17j offer a quite reasonable way out of also this predicament. The point is that 
[171 137j . the yukawa couplings of matter fermions to the MSSM Higgs receive large wave 
function corrections due to the circulation of heavy fields within loops on the lines entering 
the yukawa vertex. As a result the tree level yukawa couplings of the NMSGUT must be 
dressed before they can be matched with those in the SIMSSM : 

Y f = (1 + A T f ) ■ (Y f ) tree ■ (1 + A / )(l + A H ±) (5.22) 

Due to the large number of heavy fields the dressing of the Higgs fields can be rather large 
(>> 10). We already calculated [17] the dressing for the 10-plet component of the MSSM 
Higgs. However in our realistic fits we find that the other components (in particular those 
from the 210) can form a significant fraction of the MSSM Higgs. Above we showed that a 
completely independent line of argument requires that the doublet H be 50% derived from 
the 210-plet ! Thus the lengthy calculation of the wave function corrections for each of the 
six GUT doublets contributing to the MSSM doublet is necessary. Even from the partial 
calculation [TTJ one can see that the large value of the wave function dressing makes the 
GUT tree level matter fermion yukawa couplings (i.e {Iiab-, 9AB-, Jab} tree and therefore all 
the {y^ B )tree ) required to match the SIMSSM couplings at much smaller than they 
would be without these corrections! It is important to note that this reduction in SO(10) 
16-plet yukawa coupling magnitudes allows the d = 5 baryon violation rates - which have 
always been problematically large in super symmetric GUTs - to be reduced to acceptable 
levels r^g^Q < 10~ 34 yr _1 . The NMSGUT offers a novel and structural resolution of this 
longstanding problem by taking seriously the non trivial wave function renormalization of 
the light Higgs doublets of the MSSM by the huge number of heavy fields they are coupled 
to. Since it is the tree level couplings that enter the formulae for the inflaton dynamics in 
the full GUT it is easier to satisfy eqn. (|2.22p . Because of this and the relatively large value 
of M ~ M3 it should be be possible to achieve the required fine tuning once the full wave 
function dressing is computed. 

The embedding in the GUT has overturned our naive assumption that the lowest in- 
termediate scale would govern inflation. Instead it is rather the largest. While setting us 
the problem of finding solutions to the tuning condition, compatible both with an accurate 
fit of fermion masses and acceptable values of inflationary power spectrum and spectral 
index, it emphatically shows that the soft terms have little role to play in the fine tuning 
which belongs rather to the GUT and intermediate scale physics only. Thus the physics 
of SIMSSM driven inflation is in sharp contrast to the Dirac neutrino mass MSSM driven 
inflation[6j [7] and makes it clear that they lie counterpoised not only as regards the na- 
ture of neutrino mass but also as regards the nature of inflation and its regulating mass 
scale besides their degree of naturalness. Note that the quadratic dependence of corrections 
to soft susy parameters on the heavy masses as opposed to the logarithmic wave func- 
tion normalization of superpotential parameters makes the weaker fine tuning demands on 
superpotential parameters only in the SIMSSM case even more appealing. 
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Parameter 


K CLIUC 


Parameter 


Value 


XX 


(J.445o 






Xz 


0.142b 


M x 


7.08 x 10 


h 


1.06b x 10 


p p 

JI1J2 


2.59 x 10 , 4.405 x 10 


a 


o A A i n — 4 
Z.44 X 10 


A„ 


0.999999 


M 






n QQQQQQ 

u.yyyyyy 


p 
l 


/I 9/(q v 1 n~5 


A 

'-^tuning 


n QRQ 

u.yoy 


1 


KQ y 1 n 15 


iviX 


^ 9^ x 10 17 


A (M x ),m (M x ) 


-5.235 x 10 5 , 1.260 x 10 4 


^B(M X ) 


4.316 x 10 5 , -1.128 x 10 11 


M- 


-1.498 x 10 11 


Ml 


-1.448 x 10 11 


|Atf |,|A^ o | 
Mf 


50.254,63.930 




0.707 


4.86 x 10 13 




1.181 x 10 9 ,2.01 x 10 12 


\Vz\tree\ 


1.997 x 10- 4 


\V21tree\i \Vlltree\ 


4.489 x 10" 5 , 1.640 x 10~ 6 


Log 10 (h 2 /M) 


-18.706 


Vo, (fiend 


3.579 x 10 52 , 2.153 x 10 15 


N pi vo t, NcMB 


54.22,4.78 x 10~ 4 


A,/3 


8.82 x 10" 12 ,5.92 x 10~ 6 



Table 1: Illustrative example of relevant parameters from an accurate fit of the fermion spectrum 
in the NMSGUT which is compatible with inflationary scenario. All masses are in GeV. 
Xx,z are the accuracies of the fits to 18 known fermion mass/mixing parameters at Mx,z- 

In Table I we give an example of the relevant parameters from an accurate fit of the 
complete fermion spectrum in the NMSGUT : which has also been tuned to make it as 
compatible as possible with the inflationary scenario presented here. The complete details 
regarding the fit are given as Appendix I. It is apparent that the fine tuning between the 
yukawas proceeds as anticipated with 1 — Y = 1 = A n ^. The main problem lies in the fact 
that h 2 /M ~ 10~ 19 GeV is too large by six orders of magnitude. As a result the number 
of e- folds Nqmb is much smaller than required. However as explained the formulae used 
seriously underestimate the Higgs wave function corrections. Moreover the search of the 
huge parameter space has just begun. Thus we are confident that this problem can also be 
overcome and a completely realistic fit compatible with inflation achieved. 

Finally we remark that the single stage breaking of the simple group Spin(10) to the 
SM gauge group will lead to the formation of monopoles with a Kibble density uk ~ M x 
at the time of the GUT phase transition. However inflation by 50 or more e-folds occurring 
long after the epoch when the SO(10) monopoles are formed will dilute the monopoles to 
completely levels removing any monopole problem or signal. 

6 Discussion 

In this paper we have shown how Supersymmetric Type-I seesaw models with the typical 
superpotential couplings found in MSLRMs and MSGUTs allow an attractive and natural 
implementation of renormalizable inflection point inflation. Inflation parameters are tied to 
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seesaw parameter values and the required fine tuning is less severe and more stable than in 
the Dirac neutrino case since it is essentially independent of the supersymmetry breaking 
parameters and is governed by the physics of intermediate scales ~ 10 6 — 10 12 GeV. In the 
Dirac neutrino case [B] the opposite it is true and the inflation occurs at low scales. 

The post-inflationary reheating behaviour in the our model differs from the Dirac neu- 
trino case. The mechanism of "instant preheating" [11] applied to inflection point inflation 
models shows that oscillation after slow roll of a Susy flat direction inflaton [23] ensures ef- 
ficient transfer of all the inflaton energy into thermalized MSSM plasma within few Hubble 
times after the end of inflation and consequently a high reheat temperature T r h ~ 10 11 — 10 15 
GeV. Thus this type of model requires a gravitino mass larger than about 50TeV to re- 
main consistent with Nucleosynthesis. Such large Supersymmetry breaking scales are also 
required by the NMSGUT to fit all the fermion data[16j. The high reheat temperatures 
and the presence of the Higgs in the inflaton sit comfortably with the requirements of 
thermal [22] and non thermal Leptogenesis [12] . The current work therefore extends the al- 
ready wide scope of the New Minimal Supersymmetric GUT from a completely realistic 
theory compatible with the central paradigms of Beyond Standard Model(BSM) physics 
and predictive of parameters crucial to the discovery of Supersymmetry. It has been shown 
to potentially harbour a consistent Inflationary cosmogony tied to the central paradigms of 
seesaw neutrino mass and Leptogenesis. The complete calculation [27 \ of the wavefunction 
corrections to the tree level relations between SIMSSM and NMSGUT yukawa couplings 
will permit us to confirm the viability of our scenario in the NMSGUT context. 
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Appendix 



Parameter 


T/ n lap 


Field 

l OLy y-^ji 1 j 


A/fn t s 
(Units ofW w Gev) 

1 L/ 1 VVVO KJ J -Lkj VJ D(7 1 




0.4458 


All 1 41 


1323 58 




0.1426 


R[6 2 5/31 

JD [U, Z,, O/ OJ 


0.1912 


fll /If) -9 


25 8969 


C\R 2 11 


60 14 718 45 741 09 


722/ lu 


440 5316 

i. Lkj . 'J ' ) 1 \> 


Df3 2 7/31 


62 83 71 5 80 752 1 6 


fo, /10~ 2 


0.1066 


E\3 2 1/31 


30 46 93 55 30 


/in /10~ 6 


-0 8871 + 1599? 




55 303 830 73 942 57 

KJ KJ • KJ KJ KJ UUv. 1 KJ ^ i/liiiL' 1 


ftio /10~ 6 
"-12/ xu 


— Q 0533 -1- 8 76997 


Fll 1 21 


12 61 12 61 


"-13/ xu 


5 8257 -1-0 41 00?' 




44 84 664 62 


ZlWlO -5 


_4 iggs _i_ 19 6580i 


Gil 1 01 

u^x, -L, Uj 


048 38 75 

V/ • KJ^XKJ ^ KJ • KJ KJ ^ W • 1 KJ 


/,„, /If) -4 


8 S003 — 4 1 2557 




755 32 31 32 44 


/ioo /in- 3 

"'33/ xu 


-3 7417 - 1 65957 

KJ • 1 i -1- 1 J- • KJ KJ ij KJ V 


Ml 2 11 


1 176 40 90 62 13 


On/ID -4 

yi2/ j-« 


_n 1179 + o 09407 




1120 35 1178 92 


i/13/ iu 


5 991 1 _ n 50957 


H3 1 10/31 

J [O, -L, 1U/ OJ 


67 

KJ ,\J t 


O9./10~ 4 


9.1454 - 9.2551? 


J [3 1 4/31 


0.737, 29.75, 29.75 


A/irr 2 


-0.2982 - 0.3350? 




86.85, 798.73 


'/ 


-10.1628 + 3.9777? 


1H3 1 8/31 


100.76,972.38 


r 


0.4475 - 2.1204? 


Lf6 1 2/31 

J-/|0, ^ 5 ^/ OJ 


48 91 1571 17 

A KJ • %J J- ^ -L ^/ 1 _L ■ A 1 




0247 — 0765? 


Mf6 1 8/31 


1590.77 


s 


1.2522 + 0.4940? 


iV[6, 1,4/31 


1582.18 


s 


0.8170 + 0.8221? 


OFl, 3, 21 

L ' J 


3043.86 


m/10 16 GeV 


0.02 


P[3,3,2/3l 


21.71,2384.12 


m^/lQ m GeV 


—41 889e _iAr3 ( A ) 


G\H 3 01 


0.559 


7 


3 78 

kj- 1 O 


/?18 i oi 

JX[0, ±, UJ 


21 82 


7 


-3.5398 


S[l,3,0] 


0.9277 


X 


0.9382 + 0.6473?' 


t[3, 1,2/3] 


0.60,38.05,94.02,181.09 


Ax 


1.52 




555.93,755.27,15333.10 


A G 


-7.505 


C/[3,3,4/3] 


0.786 


Aa 3 (M z ) 


-0.004 


y[i,2,3] 


0.549 


{M^ /10 iZ GeF} 


0.001181,2.01,48.59 


W[6, 3, 2/3] 


1877.78 


{M/ 7 /10- 12 eF} 


0.3880,660.09,15968.49 


X[3,2,5/3] 


0.185,59.281,59.281 


M v (meV) 


2.148903,7.32,40.17 


^[6,2,1/3] 


0.23 


{Ev a is[f]}/irr 6 


0.025897, 44.05, 1065.71 


Z[8, 1,2] 


0.81 


Soft parameters 
at M x 

Max(\L A BCD\, \Rabcd\) 


mi = -1062.672 

fi = 4.3160 x 10 5 
M| = -1.4978 x 10 11 

9.6779 x 10~ 23 GeV -1 


m = 12603.819 

B = -1.1281 x 10 11 
M 2 = -1.4480 x 10 11 


A = -5.2347 x 10 5 

tan/3 = 50.0000 
Rbr = 1.4504 



Table 2: Fit : Column 1 contains values of the NMSGUT-SUGRY-NUHM parameters at M x de- 
rived from an accurate fit to all 18 fermion data and compatible with RG constraints. 
Unihcaton parameters and mass spectrum of superheavy and superlight fields are also 
given. The values of fi(Mx), B(Mx) are determined by RG evolution from Mz to Mx of 
the values determined by the EWRSB conditions. 
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Parameter 


1 arget = 


Uncert. = Oi 


Achieved = Oi 


t>ull = [Ui - Ui)/di 




no 1 coo 

2.031523 


0.776042 


2.036884 


0.006908 


/i a— 3 


0.990278 


0.163396 


n nor /I ot 

0.985427 


—0.029685 




0.350699 


f\ A1 A AOO 

0.014028 


n rno tct 

0.350375 


—0.023035 


/i n— 5 


7.314770 


4.264511 


8.249786 


r\ 01 00 cr r 

0.219255 


/I n — 3 


1 r 71 1 

1.385711 


0.654056 


1 r cr /i n 

1.335549 


—0.076693 


Vb 


0.438505 


0.227584 


0.496225 


0.253620 


/i n— 4 


1.190847 


0.178627 


1 100 00 

1.182832 


n A A 0^71 

—0.044871 




O A A A C A C\ 

Z.44454U 


0.366681 


2.408165 


—0.099201 


Vt 


0.519320 


0.098671 


0.529217 


0.100302 


sin 6*^2 


0.2210 


r\ AA1 r* r\r\ 

0.001600 


0.2210 


—0.0066 


sm (7^3/ 1U 


00 /loon 
zy.4zyy 


O.UUUUUU 


90 ^1 A9 

zy.oiuz 


A AI fil 
U.UlDl 


smt'23/lu 


9/1 £0^70 

34. DZ / z 


1 QAAAAA 

1.300000 


34.D44U 


A AI OA 

0.0129 




£A aoi 1 
OO.Ozll 


1 A AAAAAA 

14.U0U0U0 


KCl A/1 Q 1 

5y.y43i 


— U.U05D 


(m 12 )/lU (el/J 


/I 007Q 

4. ay < 


A ^1 OI AO 

u.oiyiuy 


4.oy (y 


a Am 
U.UU1Z 


(m 2 3J/lU (el/ J 


1.0D13 


0.31ZZ IU 


l.OOUU 


A A A A Q 

— U.UU43 


sm (7-^2 


u.zyoy 


U.UDo / oU 


A on A A 

u.zy44 


A AAA/1 

u.uuy4 


sm fl 23 


U.4olo 


U.13o5oz 


a a c;n7 
U.45y 1 


A AI 

— U.U151 


sm (7-^3 


A noco 
U.UZOZ 


n ai nnnn 
u.uiyuuu 


A A99C; 
U.UZZO 


A 1 /I QQ 

— u. i4oy 


Eigenvalues(/\u) 


0.06601 < 


0.066029 


A nc^n A A 
U.0b6U44 




Eigenvalues(A^) 


0.063539 


0.063551 


0.063566 




Eigenvalues(Ap) 


0.0 1 303 / 


A A^QA/IA 

0.0/3049 


0.0 ( 3064 




EigeTivalues(/\g) 


A AQA/1 TO 

0.0804/2 


A AQA/1 Q/1 

0.080484 


A AQA/1 AO 

0.080499 




Eigenvalues (Aq) 


A A£; icin 

0.061610 


A A£J 

0.06162z 


0.061635 




Eigenvalues^AjJ) 


A A7QKQ« 
U.U ( 0D0D 


A A79CAQ 

u.u ( ooyy 


A A70C1 1 

U.U ( oDll 




A - A „ 


Do.yoUloD 


0U.Z044 / 1 






cti 


U.04UZ + U.UUUU2 


Ql\ 


n 79on n nnnn// 
U. ( ZZU — U.UUUU? 




ct2 


n dm s j- n 091 7i 


6t2 


U.UOOI t U.U04Ut 




a 3 


-0.0405 - 0.0412i 


«3 


-0.0619 - 0.0274i 




04 


-0.6968 + 0.1200i 


CX4 


0.6213 - 0.016H 




a 5 


0.1061 + 0.0735i 


«5 


0.0585 + 0.0173i 




a 6 


0.1356 - 0.2204i 


«6 


0.1646 - 0.2294i 






0.640,0.056 


«1 , «2 


0.722,0.066 




a 3 , a 4 


0.058,0.707 


«3 , «4 


0.068,0.622 




"5 1, «6 


0.129,0.259 


«5 1, «6 


0.061,0.282 





Table 3: Fit with X x = VEjli( ^ ~ O l ) 2 /Sf = 0.4458. Target values, at M x of the fermion 
yukawa couplings and mixing parameters, together with the estimated uncertainties, 
achieved values and pulls. The eigenvalues of the wavefunction renormalization incre- 
ment matrices Aj for fermion lines and the factors for Higgs lines are given, assuming the 
external Higgs is 10-plet dominated. The Higgs fractions ai,di which control the MSSM 
fermion yukawa couplings are also given. Right handed neutrino threshold effects have 
been ignored. We have truncated numbers for display although all calculations are done 
at double precision. 
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Parameter 


SM(M Z ) 


m GUT (M z ) 


m MSSM = ( m + Am) GUT (M z ) 




2.90000 


0.75028 


3.17037 


m s /W- ? > 


55.00000 


12.14614 


51.36207 


m b 


2.90000 


3.19379 


3.05311 


m e /io~ 3 


0.48657 


0.45539 


0.46641 




0.10272 


0.09267 


0.09776 


m T 


1.74624 


1.71722 


1.71666 


m u /io~ 3 


1.27000 


1.10756 


1.28852 




0.61900 


0.53583 


0.62338 


m t 


172.50000 


143.44513 


171.73898 



Table 4: Values of standard model fermion masses in GeV at Mz compared with the masses obtained 
from values of GUT derived yukawa couplings run down from M x to Mz both before 

and after threshold corrections. Fit with X z = \Jj2Li( m f SSM ~ mf M ) 2 /(mf SSM ) 2 = 
0.1408. 



Parameter 


Value 


Parameter 


Value 


Mi 


276.93 


M Sl 


15629.87 


M 2 


942.04 


M &2 


15625.51 


M 3 


662.69 


Mu 3 

A m 
A ii 


76919.36 


M k 


3892.50 


-323296.22 


M ! 2 


283.29 




-322948.67 


M t 3 


65951.57 


^33 


-204476.65 


M Ll 


23375.27 


A 0(u) 


-391038.17 


M U 


23214.39 


A 0(u) 


-391035.62 


M U 
M-- 

di 


52427.68 
3610.56 


4 0(«) 
^33 


-211710.78 
-322645.16 


M-- 

d2 


3604.89 


40(d) 
22 


-322642.31 


M k 


134282.49 


^33 


-125043.45 




17575.24 
17572.75 


tan/3 
H{M Z ) 


50.00 
351033.09 




109825.82 


B(M Z ) 


2.4726 x 10 10 


M 2 - 


-1.1964 x 10 11 


Ml 


-1.3584 x 10 11 



Table 5: Values (GeV) in of the soft Susy parameters at Mz (evolved from the soft SUGRY- 
NUHM parameters at Mx). The values of soft Susy parameters at Mz determine the 
Susy threshold corrections to the fermion yukawas. The matching of run down fermion 
yukawas in the MSSM to the SM parameters determines soft SUGRY parameters at Mx- 
Note the heavier third sgeneration. The values of u{Mz) and the corresponding soft 
susy parameter B{Mz) = m\sm2(3/2 are determined by imposing electroweak symmetry 
breaking conditions, wia is the mass of the CP odd scalar in the in the Doublet Higgs. 
The sign of a is assumed positive. 
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Field 


Mass{GeV) 


Mg 


662.69 


M x± 


942.04,351033.11 


M x o 


276.93, 942.04, 351033.10, 351033.11 


My 


23375.180, 23214.295, 52427.637 


M- e 


3892.76, 23375.33, 277.93, 23214.55, 52422.21, 65955.95 


Mu 


15629.83, 17575.16, 15625.45, 17572.68, 76918.45, 109826.62 


M S 


3610.66, 17575.35, 3604.97, 17572.86, 109823.42, 134284.46 


M A 


1112118.78 


M H ± 


1112118.78 


M H o 


1112118.78 


M h o 


122.98 



Table 6: Spectra of supersymmetric partners calculated ignoring generation mixing effects. In- 
clusion of such effects changes the spectra only marginally. Due to large values of 
H » Mz,M\y the LSP and light chargino are essentially pure Bino and Wino(W±). 
The light gauginos and light Higgs h°, are accompanied by a light smuon and sometimes 
selectron. The rest of the sfermions have multi-TeV masses. The mini-split supersym- 
metry spectrum and large /j,, Aq parameters help avoid problems with Flavor Changing 
Neutral Currents and Charge and Color breaking/Unbounded from below(CCB/UFB) 
instability 17 . The sfermion masses are ordered by generation not magnituide. This is 
useful in identifying the spectrum calculated including generation mixing effects. Note the 
very light(right) smuon. 



Field 


Mass{GeV) 


M G 


663.15 


M x ± 


942.22,351025.61 


M x o 


276.99, 942.22, 351025.60, 351025.60 


M D 


23214.64, 23375.50, 52426.007 


Mg 


249.75, 3890.86, 23214.90, 23375.64, 52420.65, 65953.06 


Mu 


15626.70, 15631.41, 17574.07, 17576.34, 76909.50, 109817.78 


Mg 


3604.73, 3615.05, 17574.26, 17576.53, 109815.13, 134273.86 


M A 


1112398.16 


M H ± 


1112398.16 


M H0 


1112398.15 


M h o 


122.99 



Table 7: Spectra of supersymmetric partners calculated including generation mixing effects. In- 
clusion of such effects changes the spectra only marginally. Due to large values of 
(j, >> Mz,Mw the LSP and light chargino are essentially pure Bino and Wmo(W±). 
Note that the ordering of the eigenvalues in this table follows their magnitudes, compari- 
son with the previous table is necessary to identify the sfermions 
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